Exact results and new insights for models defined over small-world networks. 
First and second order phase transitions. I: General result 
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We present, as a very general method, an effective field theory to analyze models defined over 
small-world networks. Even if the exactness of the method is limited to the paramagnetic regions 
and to some special limits, it gives the exact critical behavior and the exact critical surfaces and 
percolation thresholds, and provide a clear and immediate (also in terms of calculation) insight of 
the physics. The underlying structure of the non random part of the model, i.e., the set of spins 
staying in a given lattice Co of dimension do and interacting through a fixed coupling Jq, is exactly 
taken into account. When Jo > 0, the small-world effect gives rise to the known fact that a second 
order phase transition takes place, independently of the dimension do and of the added random 
connectivity c. However, when Jo < 0, a completely different scenario emerges where, besides a spin 
glass transition, multiple first- and second-order phase transitions may take place. 

PACS numbers: 05.50.-|-q, 64.60.aq, 64.70.-p, 64.70.P- 



I. INTRODUCTION 



Since the very beginning of the pioneeristic work by 
Watts and Strogatz the interest toward small- world 
networks - an interplay between random and regular net- 
works - has been growing "exponentially" . Mainly, there 
are two reasons that have caused such a "diffusion" . 

The first reason is due to the topological proper- 
ties of the small-world network. In synthesis, if N is 
the size of the system, for any finite probability p of 
rewiring, or for any finite added random connectivity c 
(the two situations correspond to two slightly different 
procedures for building a small-world network) one has: 
a "short-distance behavior" , implying that the shortest 
distance between two arbitrarily chosen sites grows as 
1{N) ~ log(A^), as in random networks, and a large clus- 
tering coefficient, C{N) ~ 0(1), as in regular lattices. 
The interplay between these two features makes small- 
world networks representative of many realistic situations 
ranging from social networks, communications networks, 
chemical reactions networks, protein networks, neuronal 
networks, etc. 

The second reason is due the fact that, in models de- 
fined over small-world networks, despite the presence of 
an underlying finite dimensional structure - a lattice £o 
of dimension dp < oo - the existence of the short-cut 
bonds makes such models mean-field-like and - hopefully 
- exactly solvable. However, even if such a claim sounds 
intuitively correct, the complexity of these models turns 
out to be in general quite high and, compared to numeri- 
cal works, analytical results, and especially exact results, 
on small-world networks are still few i-i (for perco- 
lation an d sy nchronization problem we cite in particular 
[3 and [nf). 

In particular, for do > 1 an analytical approach seems 
to be impossible, even though, again, a mean-field critical 
behavior is expected and has been already confirmed by 



Monte Carlo (MC) simulations [l^. A natural question 
is, are we able to prove analytically such an assertion? If 
for example do = 2, does the mean- field critical behavior 
hold for any situation? Yet, does the correlation length 
diverge at the critical temperature? 

Furthermore, even if for do = 1 an exact analytical 
treatment has been performed at the level of replica sym- 
metry (RS) and one step replica symmetry breaking 
(IRSB) 0, the calculations involved are quite long and 
the solutions of the coupled equations for evaluating the 
order parameters require a certain numerical work be- 
coming rapidly hard in the IRSB case. In any case, even 
if these methods are able to give in principle exact results 
at any temperature, they are not in general suitable for 
giving a clear simple and immediate physical picture of 
the model, even - possibly - within some approximation. 
The main problem in fact stays in the presence of the 
short-loops: as soon as do > 1 these loops cannot be 
neglected and the "traditional" cavity and replica meth- 
ods seem hardly applicable. In particular, what happens, 
for example, if we set Jo negative? Shall we still expect 
a second order phase transition? And what about the 
phase diagram? 

In this paper, we present a general method to study 
random Ising models defined on small-world graphs built 
up by adding a random connectivity c over an underly- 
ing arbitrary lattice having dimension dp. We shall 
then show that this method provides - in a very simple 
and physically appealing way - the answers to the above 
questions and many others. 

As an effective field theory, roughly speaking, the 
method generalizes the mean-field equation m = 
tanh(/3Jm) to take into account the presence of the short- 
range couplings Jo besides the long-range ones J. As we 
will show, the magnetization m of the model defined over 
the small-world network, shortly the random model, be- 
haves as the magnetization mo of the model defined over 
Cq, shortly the unperturbed model, but immersed in an 
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effective field to be determined self-consistently. Even if 
the exactness of this method is limited to the paramag- 
netic regions (P), it gives the exact critical behavior and 
the exact critical surfaces, and provide simple qualitative 
good estimations of the correlation functions in the ferro- 
magnetic (F) and spin glass regions (SG). Furthermore, 
in unfrustrated systems, the method becomes exact at 
any temperature in the two limits c — > 0^ and c — > oo. 

The consequences of such a general result are remark- 
able from both the theoretical and the practical point of 
view. Once the explicit form of the magnetization of the 
unperturbed model, mo = toq (/3 Jo , , as a function of 
the couplings Jo and of the external field h is known, an- 
alytically or numerically, we get an approximation to the 
full solution of the random model analytically or numer- 
ically, respectively, becoming exact in the P region. If 
we do not have mo — mo{PJo, Ph) but we know at least 
some of its properties, we can still use these properties 
to derive certain exact relations and the critical behavior 
of the random model. 

In this paper (part I), after discussing the self- 
consistent equations, we will mainly focus on applying 
them to study the critical surfaces and the critical be- 
havior in general, whereas in a forthcoming paper (part 
II), after showing some examples, we will apply the 
method to study several models of interest which can 
be solved analytically (and very easily) as for them we 
know mo(/3Jo, Ph). 

We stress that, the critical surfaces, as well as the 
correlation functions in the P region provided by our 
method, are exact and - quite interestingly - not based 
on any special ansatz as the replica-symmetry and the 
tree-like ansatz. We prove in particular that: for Jo > 0, 
independently of the added random connectivity c and 
of the underlying dimension Joi we always have a second 
order phase transition with the classical mean-field crit- 
ical indices but with a finite correlation length if calcu- 
lated along the distance defined by the underlying lattice 
£0; whereas for Jo < 0, we show that, as soon as c is 
sufficiently large, there exist at least two critical temper- 
atures which, depending on the behavior of xo(/3Jo, Ph) - 
the susceptibility of the unperturbed system - correspond 
to first or second order phase transitions. 

The paper is organized as follows. In Sec. II we in- 
troduce the class of small-world networks over which we 
define the random Ising models, stressing some impor- 
tant difference concerning the definition of the correla- 
tion functions with respect to the definition of the corre- 
lation functions one usually considers in "ordinary" ran- 
dom models. In Sec. Ill we present our method: in Sec. 
IIIA we provide the self-consistent equations and their re- 
lations with physical correlation functions, in Sec. IIIB 
we analyze the stability of the solutions of the self con- 
sistent equations and the critical surface and behavior of 
the system. We separate the Sec. IIIB in the sub-cases 
Jo > and Jo < 0. In Sec. IIIC we discuss the limits 
of the method. In Sec. HID we study the stability be- 
tween the F and the SG phases and the phase diagram. 



Finally, in Sec. HIE we mention how to generalize the 
method to cases in which there are more different short- 
range couplings Jo, and how to generalize the method 
to analyze possible disordered antiferromagnetism. The 
successive Sees. IV, V and VI are devoted to the deriva- 
tion of the method. The starting point of the proof is 
given in Sec. IV and is based on a general mapping be- 
tween a random model and a non random one tlSj-^lsl 
suitably adapted to our case. The self-consistent equa- 
tions are then easily derived in Sec. V. Note that, apart 
from the equations concerning the stability between the 
P-F and the P-SG transitions, which are derived in Sec. 
VI, the derivations of the equations provided in the Sec. 
IIIB are mostly left to the reader, since they can be easily 
obtained by standard arguments of statistical mechanics 
using the Landau free energy ■0(m) that we provide and 
that is derived in Sec. V too. Finally, in Sec. VII we 
draw some conclusions. In Appendix A we generalize the 
method to inhomogeneous external fields to make clear 
the subtle behavior of the correlation functions in small- 
world models. 



II. RANDOM ISING MODELS ON 
SMALL- WORLD NETWORKS 

The family of models we shall consider are random 
Ising models constructed by super-imposing random 
graphs with finite average connectivity c onto some given 
lattice £0 whose set of bonds and dimension will 

be indicated by Fo and Jo, respectively. Given an Ising 
model {the unperturbed model) of N spins coupled over 
Cq through a coupling Jo with Hamiltonian 

Ho = -Jo ^ a^aj-h^ai, (1) 

(ij)ero * 

and given an ensemble C of unconstrained random graphs 
c, c G C, whose bonds are determined by the adjacency 
matrix elements Cij = 0, 1, we define the corresponding 
small-world model as described by the following Hamil- 
tonian 

Hc,J = Hq-'^ Cy JijCTidj , (2) 



the free energy F and the averages (C)' being defined in 
the usual (quenched) way as (/? — 1/T) 

-PF = Y, P{c) I dV ({J,,,}) log (Ze,j) (3) 

and 

where Z(.^j is the partition function of the quenched 
system 
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(5) 



{0)c,j the Boltzmann-average of the quenched system 
((O) depends on the given reahzation of the J's and of 
c: (O) — {0)c;j', for shortness we later will omit to write 
these dependencies) 



(0>c,. 



Z. 



(6) 



and dV {{Jij}) and P{c) are two product measures 
given in terms of two normalized measures diJ,{Jij) > 
and p{cij) > 0, respectively: 

dP({J,.,})= Y[ dAi(J„), yd/x(J,,,) = l, (7) 



P{c)= \[ p{c,,,), P(c*,j) = l- 



(8) 



The variables Cij € {0, 1} specify whether a "long- 
range" bond between the sites i and j is present j = 1) 
or absent {cij = 0), whereas the Jij 's are the random 
variables of the given bond (i, j ) . For the j 's we will not 
assume any particular distribution, while, to be specific, 
for the Cij 's we shall consider the following distribution 



N 



(9) 



This choice leads in the thermodynamic limit N ^ oo to 
a number of long range connections per site distributed 
according to a Poisson law with mean c > (so that in 
average there are in total cN/2 bonds). Note however 
that the main results we report in the next section are 
easily generalizable to any case in which Eq. ([5]) holds, or 
holds only in the thermodynamic limit due a sufficiently 
small number of constrains among the matrix elements 

The quantities of major interest are the averages, and 
the quadratic averages, of the correlation functions which 
for shortness will be indicated by C^^^ and C*^^^ . For ex- 
ample, the following are non connected correlation func- 
tions of order k: 



C(2) 



(10) 

(11) 



where fc > 1 and the indices ii, . . . ,ik are supposed all dif- 
ferent. For shortness we will keep on to use the symbols 
C*^^-* and C^^-* also for the connected correlation function 
since, as we shall see in the next section, they obey to 
the same rules of transformations. We point out that 
the set of indices ii, . . . ,ik is fixed along the process of 



the two averages. This implies in particular that, if we 
consider the spin with index i and the spin with index j, 
their distance remains undefined, or more precisely, the 
only meaningful distance between i and j, is the distance 
defined over Cq, i.e., the Euclidean distance between i 
and j, which we will indicate as ||i ~ j H^j. 

Therefore, throughout this paper, it must be kept in 
mind that, for example, C^^dji — = {(Ticrj) is very 
different from the correlation function G^^^{1) of two 
points at a fixed distance /, / being here the distance 
defined over both Cq and the random graph c, i.e., the 
minimum number of bonds to join two points among both 
the bonds of Fq and the bonds of the random graph c. 
In fact, if, for Jq = 0, one considers all the possible re- 
alizations of the Poisson graph, and then all the possible 
distances I between two given points i and j, one has 



N 



= ^p^(OgW(0 



(12) 



1=1 



where here P;v(0 the probability that, in the system 
with N spins, the shortest path between the vertices i 
and j has length I. 

If we now use G'-^\l) - (tanh(/3J))' 16] (in the P 
region holds the equality) and the fact that the average 
of I with respect to Pn{1) is of the order log(iV), we see 
that the two point connected correlation function 
goes to in the thermodynamic limit. Similarly, all the 
connected correlation functions defined in this way are 
zero in this limit. Note however, that this independence 
of the variables holds only if Jq = 0. This discussion will 
be more deeply analyzed along the proof by using another 
point of view, based on mapping the random model to a 
suitable fully connected model. 



III. AN EFFECTIVE FIELD THEORY 

A. The self-consistent equations 

Depending on the temperature T, and on the param- 
eters of the probability distributions, dfi and p, the ran- 
dom model may stably stay either in the P, in the F, or in 
the SG phase. In our approach for the F and SG phases 
there are two natural order parameters that will be indi- 
cated by TO^^-* and m^^'^\ Similarly, for any correlation 
function, quadratic or not, there are two natural quanti- 
ties indicated by C^^' and C*^^*^', and that in turn will 
be calculated in terms of m*-^) and m^^'^\ respectively. 
To avoid confusion, it should be kept in mind that in our 
approach, for any observable O, there are - in principle 
- always two solutions that we label as F and SG, but, 
as we shall discuss in Sec. HID, for any temperature, 
only one of the two solutions is stable and useful in the 
thermodynamic limit. 
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In the following, we will use the label o to specify that 
we are referring to the unperturbed model with Hamilto- 
nian ([1]). Let mo{pJo, (3h) be the stable magnetization of 
the unperturbed model with coupling Jo and in the pres- 
ence of a uniform external field h at inverse temperature 
p. Then, the order parameters rrS^\ I]=F,SG, satisfy 
the following self-consistent decoupled equations 

m(^) = maip/p ,l3J^^'^m'^^^ + j3h), (13) 

where the effective couplings J^^\ J^^'^\ Jq^^ and Jg^*^^ 
are given by 



dfj,{Jij) tanh(/3Ji^ 



(14) 



S=SG, the true solution m^'^\ or leading solution, is the 
one that minimizes L^^': 



min L^^) (m) . (21) 

mG[-l,l] 



Finally, let k be the order of a given correlation func- 
tion C(i) or C(2). The averages and the quadratic aver- 
ages over the disorder, C*^^^ and C^^\ are related to C^^-* 
and C*^^*^-* , as follows 

C(i) = C<^\ inF, (22) 

C(i) = 0, kodd, inSG, (23) 

^ fjiSG)^ keven, in SO, (24) 

and 



/3J(SG) d^(Jij)tanh2(/3J,j), 



7(F) 



and 



P4^^^ = tanh-i(tanh2(/3Jo)). 



(15) 



(16) 



(17) 



Note that '\ > Jq 



(F)| .(SG) 

For the correlation functions C'-^\ E=F,SG, for suffi- 
ciently large N we have 



C(^) = Co(/3jf \/3J(^)m(^) +0 (^1) , (18) 

where Cq {PJo i Ph) is the correlation function of the un- 
perturbed (non random) model. 

Concerning the free energy density / we have 

= -| j d^x{.h^,)\og[coMfiJ^,,)] 

hm i- V log[cosh(/3Jo)]-log[2cosh(/3/i)] 

/ — >on /V ' * 



N-*oo N 

1 



(ij)ero 



{J™Ll^ E l°g cosh(/34 



(ij)6ro 

-f log[2cosh(/3/i)]} X y + iL(^)(m(^)), 
where I = 1,2 for E=F,SG, respectively, and 



(19) 



L(^)(r 



/?/o(/3#\/5j(''^m + /3/i)(^0) 



fo{(3Jo, f3h) being the free energy density in the thermo- 
dynamic limit of the unperturbed model with coupling 
Jo and in the presence of an external field h, at inverse 
temperature /3. 

For given /3, among all the possible solutions of Eqs. 
P3p . in the thermodynamic limit, for both S=F and 



= , i 

^(2) ^ (^ciSG)y 



in F. 



in SG. 



(25) 
(26) 



From Eqs. (^5]) and (|^ for fc = 1, we note that the 
Edward- Anderson order parameter C^^-* = (u)-^ = (/ba 
is equal to (C^^^))^ = (to^^^))^ only in the SG phase, 
whereas in the F phase we have qea = (m^F))^, There- 



fore, since to^^G-* ^ m'-^-', m^^*^' is not equal to y/qEA] 
in our approach m^^*^) represents a sort of a spin glass 
order parameter. 

The localization and the reciprocal stability between 
the F and SG phases will be discussed in Sec. HID. 
Note however that, at least for lattices Cq having only 
loops of even length, the stable P region is always that 
corresponding to a P-F phase diagram, so that in the P 
region the correlation functions must be calculated only 
through Eqs. ^ and 

As an immediate consequence of Eq. (IT3l) we get the 
susceptibility x'^-* of the random model: 



X 



1 



, (27) 



where xo stands for the susceptibility xo of unper- 
turbed model divided by f3 (we will adopt throughout 
this dimensionless definition of the susceptibility): 



dmo (/3Jo, Ph) 1 dmo {(3 Jo, I3h) 



d{ph) 



dh 



(28) 



and similarly for the random model. 

For the case S =F without disorder {dfi{J') = 5 {J' — 
J)dJ'), Eq. (|27|) was already derived in p] by series 
expansion techniques at zero field (h = 0) in the P region 
(where m = 0). 

Another remarkable consequence of our theory comes 
from Eq. (fT5|) . We see in fact that in the thermodynamic 
limit any correlation function of the random model fits 
with the correlation function of the unperturbed model 
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but immersed in an effective field that is exactly zero in 
the P region and zero external field {h = 0). In other 
words, in terms of correlation functions, in the P region, 
the random model and the unperturbed model are in- 

(SG) 

distinguishable (modulo the transformation Jq ^ Jg 
for S =F). Note however that this assertion holds only 
for a given correlation function calculated in the ther- 
modynamic limit. In fact, the corrective 0(1/-/V) term 
appearing in the rhs of Eq. (fTS]) cannot be neglected 
when we sum the correlation functions over all the sites 
i € Co, as to calculate the susceptibility; yet it is just this 
corrective 0(1/A^) term that gives rise to the singularities 
in the random model. 

More precisely, for the two point connected correlation 
function: 



where Z = 1, 2 for E = F, SG, respectively, if 



(29) 



(30) 



we have 

= Xo;^,,(/3jf\/3j(^^™(^^+/3M + 



N 



(31) 



Eq. ([31]) clarifies the structure of the correlation func- 
tions in small-world models. In the rhs we have two 
terms: the former is a distance-dependent short-range 
term whose finite correlation length, for T ^ T^q (T^q 
being the critical temperature of the unperturbed model 
with coupling Jq ), makes it normalizable, the latter 
is instead a distance-independent long-range term which 
turns out to be normalizable thanks to the 1/A^ factor. 
Once summed, both the terms give a finite contribution 
to the susceptibility. It is immediate to verify that by 
summing xl / over all the indices i,j G Cq and dividing 
by N we get back - as it must be - Eq. Eq. ([51]) 

will be derived in Appendix A where we generalize the 
theory to a non homogeneous external field. 



B. Stability: critical surfaces and critical behavior 



Note that, for /3 sufficiently small (see later), Eq. (fT3|) 
has always the solution m^^-* = 0, and furthermore, if 
m*^^^ is a solution, — m^^^ is a solution as well. From 
now on, if not explicitly said, we will refer only to the 
positive (possibly zero) solution, the negative one being 
understood. A solution m'^^ of Eq. (|13p is stable (but 
in general not unique) if 



Xo 



)m(^) +(3h) > 0. (32) 



For what follows, we need to rewrite the non trivial 
part of the free energy density L^^)(m) as 

L(^)(m) - /3/o(/3jf\o) -TOo(/3jf\o)/3/i 

+^^^'^ (m) , (33) 

where the introduced term plays the role of a Lan- 
dau free energy density and is responsible for the critical 
behavior of the system. Around m = 0, up to terms 
0{h?) and Oim^h), ij^'^^{m) can be expanded as follows 



where 



-fA(/3/o)(/3jf\/3j(^)m), (34) 
-/3./(^^Xo(/?jf\0)]/3j(^), (35) 



d{l3h)- 



d{(3hy 



■Xo 



■Xo 



[(3jP,l3h 



3! 



5! 



(36) 



(37) 



h(^) ^ mo \0) J(^) -f xo \0) J^^^h, (38) 
finally, the last term A (/?/o) (/3Jo, Z?^^^^ 

m) is defined 

implicitly to render Eqs. and ([M]) exact, but terms 
0{h'^) and 0{m^h), explicitly: 



A 

CK 

E 



Q2k-2 

d{J3h}^ 



2fe-2 



XO 



t3h=0 



-<39) 



We recall that the /c — 2-th derivative of xo (^PJo^'^ , Ph 
with respect to the second argument, calculated at 
h — 0, gives the total sum of all the k-ih cumu- 

lants normalized to N: dpJ^'^xo{pJo^\P^\h=o = 

{^ii ■ ■ ■ '^iJo^/N, where {a,, ■ ■ ■ a,^)'^^^ stands 
for the cumulant, or connected correlation function, of or- 
der k of the unperturbed model, {(Ji^ai^)'"f^^ = {oi^ai^)o — 
(o'ii)o(o'i2)o, etc.. Note that, apart from the sign, these 
terms are proportional to the Binder cumulants [iTj 
(which are all zero above Tco for fc > 2) only for N finite. 
In the thermodynamic limit the terms h^'^\ c'-^\ . . ., in 
general are non zero and take into account the large de- 
viations of the block-spin distribution functions from the 
gaussian distribution. 
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Let Tc ' = 1/Pc be the critical temperatures, if any, 
of the random model and let t^^^^ be the corresponding 
reduced temperatures: 



t 



T -T, 



/3f ^ - /? 



0(i 



(40) 



Here, the term "critical temperature" , stands for any 
temperature where some singularity shows up. However, 
if we limit ourself to consider only the critical tempera- 
tures crossing which the system passes from a P region 
to a non P region, from Eq. ([5^ it is easy to see that, 
independently on the sign of Jo and on the nature of the 
phase transition, we have the important inequalities 



(41) 



where we have introduced P^q , the inverse critical tem- 

perature of the unperturbed model with coupling Jq 
and zero external field. If more than one critical tem- 
perature is present in the unperturbed model, is the 
value corresponding to the smallest value of these crit- 
ical temperatures (highest in terms of /?). Formally we 
set /3^Q — oo if no phase transition is present in the un- 
perturbed model. A consequence of Eq. ([55)1 is that, in 
studying the critical behavior of the system for h = 0, 
we can put /i'^^ = 0. Throughout this paper, we shall 
reserve the name critical temperature of the unperturbed 
model as a P-F critical temperature through which the 
magnetization mo(/3Jo:0) passes from a zero to a non 
zero value, continuously or not. This implies, in particu- 
lar, that for Jo < we have - formally - (3cQ — oo. 

In this paper we shall study only the order parameters 
m^^-* and m^^'~^-' , whereas we will give only few remarks on 
how to generalize the method for possible antiferromag- 
netic order parameters. We point out however that the 
existence of possible antiferromagnetic transitions of the 
unperturbed model does not affect the results we present 
in this paper. 

It is convenient to distinguish the cases Jo > and 
Jo < 0, since they give rise to two strictly different sce- 
narios. 



1. The case Jo > 

In this case /JJ^^^Xo {l3Jo,l3J'-^^m''^'^ + I3h) is an in- 

creasing function of /? for /3 < /3^q (and for h> 0). As a 
consequence, we have that for sufficiently low tempera- 
tures, the solution m^^^ = of Eq. (fT5|) becomes unsta- 
ble and two - and only two - non zero solutions im^^^ are 

(F) 

instead favored. The inverse critical temperatures /3c 
and Pc^'^'^ can be determined by developing - for h — - 
Eqs. for small m^^^^ and m'^^'^\ respectively, which, 
in terms of xo gives the following exact equation 



Xo (/3f)jf\0)/3f)j(^ 



1, 



/3f ^ < 



(42) 



where the constrain /3c < P^q excludes other possible 
spurious solutions that may appear when do > 2. 

The critical behavior of the system can be derived by 
developing Eqs. for small fields. Alternatively, one 
can study the critical behavior by analyzing the Landau 
free energy density ?/j(^)(m*^^)) given by Eq. 

In the following we will suppose that for Jo > 0, 
be positive (we have checked this hypothesis in all the 
models we have until now considered and that will be 
reported in the forthcoming part II of the work). Fur- 
thermore, even if the sign of c^^^ cannot be in general 
a priori established, for the convexity of the function 
/o with respect to (3h, the sum of the six-th term with 

A(/3/o) (/3J,$^\/3j(^)m(^)), in Eq. ^ must go nec- 
essarily to +00 for TO*^^^ — > oo. In conclusion, when 
Jo > 0, for the critical behavior of the system, the only 
relevant parameters of -0'^^ are a'^\ fe*^^-' and /i^^-' = 
x{f3JQ^\o)J^^^h, so that the critical behavior can be 
immediately derived as in the Landau theory for the so 
called m** model. On noting that 

a(^' > 0, for t(^) > 0, 
a(^) < 0, for t(^) < 0, 



(43) 



it is convenient to define 



d_ 



so that we have 



(44) 



(45) 



Note that, due to the fact that Jo > 0, A^^' > 0, 
and, as already mentioned, > as well. By using 
Eq. (|45p for [3 < Pco and near (3c , we see that the 



minimum m*-^) of V'^^'^i the solution of Eq. ([Te 
near the critical temperature, is given by 

t(^) >0, 

(46) 
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I _ +o(t(^)), i(^) <0. 



Similarly, we can write general formulas for the sus- 
ceptibility and the equation of state. We have 

/3j'^)xo(/3jf\0) ^ 



.71^7+0(1), t(^)>0 



(47) 



m(^)(/i) 



+0 



(/.i)(.48) 



Finally, on using Eqs. (|34p and we get that the 
specific heat C*^^-' has the following finite jump disconti- 
nuity at /3c 

C(S) 



t(^) >0, 



(49) 
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where Cc is the continuous part of the specific heat cor- 
responding to the part of the free energy density without 

Hence, as a very general result, independently of the 
underlying dimension dp and the added random connec- 
tivity c, provided positive, we recover that the random 
model has always a mean-field critical behavior with a 
second order phase transition with the classical expo- 
nents P = 1/2, -f = -f' — 1, S — 3 and a ~ a' ^ 0, 
and certain constant coefficients depending on the sus- 

ceptibility xo and its derivatives calculated at P = (3c 
and external field h ~ 0. Note however, that the corre- 
lation length of the system calculated along the distance 

of £oj II ■ Ho J remains finite also at /3c • In fact, from Eq. 



(ITS)) , for the two point correlation function at distance 
r = ||i — illjj in £0 we have 

C(^)(r) = Co{(34^\ pJ^^^m^^^ + (3h;r). (50) 

If we now assume for Co(/3Jo, 0; r) the following general 
Ornstein-Zernike form 



CoiPJo,0;r) 



(51) 



/o('') = foiPJo'Tf) being a smooth function of r (which 
has not to be confused with the free energy density), and 
Co = ^oiPJo) the correlation length, which is supposed to 
diverge only at /3cO (if any), on comparing Eqs. ([50]) and 
(EI]) for /3 > plP^ we have 



C(^)(r) 



where 



and 



/^'^^(0 = /o(/3jf^r), 
C^^'=eo(/3jf^). 



Therefore, due to the inequalities (f4T|) . we see that 



(52) 

(53) 
(54) 

(55) 



The knowledge of Co(/3Jo, (ih] r) also for h Q would al- 
lows us to find the general expression for C*^^-* (r) through 

Eq. ([50]) also for /? > /3c ■ However, since Coi(3Jo,l3h; r) 
has no critical behavior for ft, ^ 0, it follows that C*^^^ (r) 
cannot have a critical behavior for (3 > (3c either (and 

then also for (3 ^ (3c from the right). This result is 
consistent with [2]. 



2. The case Jo < 

In this case Jq < 0, so that - in general - 
/Jj'^^Xo {l3Jo, (3J'^^^mS^^) is no longer a monotonic func- 
tion of (3. However, it is easy to see that that for /3 = 



and (3 00, this function goes to 0. Therefore, for a 
sufficiently large connectivity c, from Eq. we see 

that it may there appear at least two regions where the 
paramagnetic solution m*^^) = is stable, separated by 
a third region in which a non zero solution is instead 
stable. However the situation is even more complicated 
since, unlike the case Jq > 0, the non monotonicity of 
(3J^'^^Xo{PJo:PJ''^^m'^^^) reflects also in the fact that 
the self-consistent Eq. (fT3| for S =F may have more 
solutions of the kind ±m^^\ ±m'{F), . . . which are still 
stable with respect to the stability condition (|32p . for 
h = 0. We face in fact here the problem to compare more 
stable solutions. According to Eq. ((2T|) . in the thermo- 
dynamic limit, among all the possible stable solutions, 
only m^^\ the solution that minimizes L^^\ survives, 
whereas the not leading ones play the role of metastable 
states. This kind of scenario, which includes also finite 
jump discontinuities, has been besides observed in the 
context of small-world neural networks in [18] where we 
even observe some analogy in the used formalism, at least 
for the simplest case of one binary pattern. 

From Eqs. and ([57)1 we see that the signs of 

the Landau coefhcients a(^\ b'^^\ c(^\ . . ., are func- 
tions of /3 and Jq only. Given Jq < 0, the most im- 
portant quantity that features the non monotonicity of 
PJ^^'^Xo (/3Jo, /3J*-^-'m^^)) is the minimum value of /3 over 
which b^^'' becomes negative: 



(56) 



(57) 



6(^)<0, /3>/3f\ 
or, in terms of temperatures 

0, T<T^^\ 

The equation for /3, , as a function of Jq, defines a 
point where fc'^-* — 0. If Jq < 0, the most general equa- 
tion for a generic critical temperature is no longer given 
by Eq. (UH). In fact, in general, a critical temperature 
now is any temperature where the stable and leading so- 
lution rh^^^ may have a singular behavior, also with finite 
jumps between two non zero values. 

There are some simplification when for the Landau co- 
efficient c^^', we have c'^' > 0, or at least c*^^-' > out 
of the P region (this is the case of a model that we will 
report in the part II of the work). In this situation in 
fact, from Eq. ([M)) we see that a^^\ fo*^^-* and c*^^^ are the 
only relevant terms for the critical behavior of the system 
and - for small values of m^^^ - we can again apply the 
Landau theory, this time for the so called model. In 
such a case, for the solution ?7i(^) we have 



,(F) 



[\/{b(^^y-^ai^^c^^^-b(^^y if. 



a'^^^ < 0, or 

a'^^'> > and b'^^^ < -4 



(58) 
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whereas 



i(^) = 0, if, 



a(F) > and b'^) > (59) 



From Eqs. ^ and ^ we see that, if 6^^) > 0, we 
have a second order phase transition and Eqs. (|42p - ([55|) 
are recovered with Eq. ([55]) becoming the second of Eqs. 
(|46p for small and negative values of a'^^. However, from 
Eq. ([55)) we see that, if b'^^'> is sufficiently negative, we 
have a first order phase transition which, for small values 
of a^^\ gives 



6(F) 



f - 



a^"^) < and 6'^) < 0, or 



if, 



/„(F)„(F) 

a(F) > and 5^^) < (60) 



From Eq. (US]) we see that the line b^^) -4v^o(f)c(^')/3 
with a^^' > establishes a line of first order transitions 
over which m*-^) changes discontinuously from zero to 



Am(F) 



n(F) 



(61) 



The point a'-^^ = b'-^^ = is a tricritical point where 
the second and first order transition lines meet. If we 
approach the tricritical point along the line fc^^^ we 
get the critical indices a ~ 1/2, a' = 0, /3 = 1/4, 7 = 
7' = 1 and (5 = 5. However, this critical behavior along 
the line &(^) = has not a great practical interest since 
from Eq. ([36|l we see that it is not possible to keep 6^^) 
constant and zero as the temperature varies. Finally, 
we point out that, even if c^^) > 0, when the transition 
is of the first order, Eqs. ((58)) and (|60)) hold only for 
b^^\ and then a^^\ sufficiently small, since only in such 
a case the finite discontinuity of tti^^) is small and then 
the truncation of the Landau free energy term ip'^^'> to a 
finite order meaningful. Note that this question implies 
also that we cannot establish a simple and general rule to 
determine the critical temperature of a first order phase 
transition (we will return soon on this point). 

When c^^) < 0, the Landau theory of the model 
cannot be of course applied. However, as in the case 
Jo > 0, even if the sign of c^^-* cannot be a pri- 
ori established, for the convexity of the function /o 
with respect to (3h, the sum of the six-th term with 

A(/3/o) i^f34^\pj'^^'>m^^^y in Eq. dH must go nec- 
essarily to +00 for m*-^) — > 00 and a qualitative similar 
behavior of the model is expected. In general, when 
Jo < 0, the exact results are limited to the following ones. 

From now on, if not otherwise explicitly said, we shall 
reserve the name critical temperature, whose inverse 
value of (3 we still indicate with (3c , to any tempera- 
ture on the boundary of a P region (through which m^^) 



passes from to a non zero value, continuously or not). 
For each critical temperature, depending on the value of 
(3i , we have three possible scenario of phase transitions: 



/3f ) > /3f ) 



first order, 
tricritical point, 
second order. 



(62) 



or, in terms of temperatures, 



rpiF) ^ y(F) 

tP > rh 



first order, 
tricritical point, 
second order. 



(63) 



Note that, according to our definition of critical tem- 
perature, the critical behavior described by Eqs. (|55)) - 
(1611) represents a particular case of the general scenario 
expressed by Eqs. (|62p. We see also that, in general, for 
the critical exponent /3 we have /3 < 1/4. 

In the case in which corresponds to a second order 
phase transition, or in the case in which a'^^ < out of 
the P region (at least immediately near the critical tem- 
perature), f3c can be exactly calculated by Eq. 
When we are not in such cases, the only exact way to 
determine the critical temperature is to find the full so- 
lution for m^^) which consists in looking numerically for 
all the possible solutions of Eq. (fT3)) and - among those 
satisfying the stability condition (|5^ - selecting the one 
that gives the minimum value of L^^\ 



C. Level of accuracy of the method 



In the P region, Eqs. ()13ll3ip are exact, whereas in the 
other regions provide an effective approximation whose 
level of accuracy depends on the details of the model. In 
particular, in the absence of frustration the method be- 
comes exact at any temperature in two important limits: 
in the limit c — > 0+, in the case of second-order phase 
transitions, due to a simple continuity argument; and in 
the limit c — > cxd, due to the fact that in this case the 
system becomes a suitable fully connected model exactly 
described by the self-consistent equations (fT^ (of course, 
when c 00, to have a finite critical temperature one 
has to renormalize the average of the coupling by c). 

However, for any c > 0, off of the P region and in- 
finitely near the critical temperature, Eqs. (|13m8p are 
able to give the exact critical behavior in the sense of 
the critical indices and, in the limit of low temperatures, 
Eqs. ()13m7p provide the exact percolation threshold. In 
general, as for the SK model, which can be seen as a 
particular model with Jo = 0, the level of accuracy is 
better for the F phase rather than for the SG one and 
this is particularly true for the free energy density f^^\ 
Eq. (fTO)) . In fact, though the derivatives of /(^^ are 
expected to give a good qualitative and partly also a 
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quantitative description of the system, Z*-^*^-' itself can 
give wrong results when the SG phase at low tempera- 
tures is considered. We warn the reader that in a model 
with Jo — 0, and asymmetrical distribution dfi{Jij) with 
variance J, the method gives a ground state energy per 
site u^^'^\ which grows with c as it^^*^) ~ —Jc, whereas 
the correct result is expected to be u^^*^^ ~ —J^/c. As 
a consequence, in the SK model, in the limit /3 — s- oo, 
the method gives a completely wrong result with an infi- 
nite energy. We stress however that the order parameters 
m'^^ and m*^^^', and then also the correlation functions, 
by construction, are exact in the zero temperature limit. 



that in this region there may be also mixed phases and 
re-entrance phenomena) . 

In Sec. VI we prove that there are four possible kind of 
phase diagrams that may occur according to the cases 
(1) (Jo > 0; do < 2, or do = oo), (2) (Jo > 0; 2 < 
do < oo), (3) (Jo < 0; do < 2, or do = oo), and (4) 
(Jo<0; 2 < do < oo). The four kind of possible phase 
diagrams are schematically depicted in Figs. [l]|l] in the 
plane (T, X). 

1. Jo > 



D. Phase diagram 

The inverse critical temperature (3c of the random 
model is in general a non single- value function of X: 
Pc = Pci^), where X represents symbolically the pa- 
rameters of the probability dfj, for the couplings Ji j-, and 
the parameter c, the average connectivity (which is also a 
parameter of the probability distribution of the short-cut 
bonds). The parameters of d/x can be expressed through 
the moments of d/i, and as they vary the probability d/i 
changes. For example, if dfj, is a Gaussian distribution, 
as in the SK model, there are only two parameters given 
by the first and second moment. 

In the thermodynamic limit, only one of the two so- 
lutions with label F or SG survives, and it is the so- 
lution having minimum free energy. In principle, were 
our method exact at all temperatures, we were able to 
derive exactly all the phase diagram. However, in our 
method, the solution with label F or SG are exact only 
in their own P region, i.e., the region where m*-^) = 
or TO*-^^-* = 0, respectively. Unfortunately, according to 
what we have seen in Sec. IIIC, whereas the solution 
with label F is still a good approximation also out of 
the P region, in frustrated model (where the variance of 
d/i is large if compared to the first moment) the free en- 
ergy of the solution with label SG becomes completely 
wrong at low temperatures. In conclusion, therefore, we 
are not able to give in general the exact boundary be- 
tween the solution with label F and the solution with 
label SG, and in particular we are not able to give the 
frontier F/SG. However, within some limitations which 
we now prescribe, we are able to give the exact critical 
surface, i.e., the boundary with the P phase, establish- 
ing which one - in the thermodynamic limit - of the two 
critical boundaries, P-F or P-SG, is stable (we will use 
here the more common expression "stable" instead of the 
expression "leading" ) , and to localize some regions of the 
phase diagram for which we can say exactly whether the 
stable solution is P, F, or SG. We will prove the stabil- 
ity of these solutions in Sec. VI. When for a region we 
are not able to discriminate between the solution with 
label F and the solution with label SG and they are both 
out of their own P region, we will indicate such a region 
with the symbol "SG and/or F" (stressing in this way 



As we have seen in Sec. IIIBl, if Jo > 0, for both the 
solution with label F and SG, we have one - and only one 
- critical temperature. 

The stable inverse critical temperature (3c satisfies the 
following rules. 

Case (1): If do < 2 and Jo is a finite range coupling, 
or else do = oo at least in a broad sense (see [lal), (^c{^) 
is a single- value function of X, and we have 



/3, = min{/3f),/3f«)}, 
or, in terms of temperatures 

Te = max{Tf),rfG)}. 



(64) 



(65) 



A schematic representation of this case is given in Fig. 1. 
Case (2): If instead 2 < do < oo we have 



or, in terms of temperatures 



if /3f°)>/3r\ 
if /3f^)</5r, 



(66) 



rp _ rTi(F) 



if t^^^UtP, 
if Tf«)>rf). 



(67) 

Notice in particular that the second line of Eq. ([5^ (or 
Eq. ([57]) ) does not exclude that Pc{X) may be a non 
single- value function oi X. A schematic representation 
of this case is given in Fig. 2. 

2. Jo < 

As we have seen in Sec. IIIB2, if Jo < 0, for a suffi- 
ciently large connectivity c, the solution with label F has 
at least two separated P regions corresponding to two 
critical temperatures. Here we assume that the underly- 
ing lattice Co has only loops of even length so that, for 
example, triangular lattices are here excluded. Let us 
suppose to have for the solution with label F only two 
critical temperatures (the minimum number, if Jo < 0), 
and let be 



(68) 



10 



or, in terms of temperatures, 

^.(F) ^(F) 



(69) 



In general we have the foUowing scenario. 

Case (3): If do < 2 and Jq is a finite range coupling, 
or do = oo in a broad sense (see [l^), /?c2(-X') is a single- 
value function of X and satisfies Eq. (or, in terms 
of temperatures, Eq. for Tc2). The other critical 

inverse temperature /3ci(-X') is instead: either a two- value 
function of X and we have 



•cl 



«(SG) 



(3c 



if 



(70) 



(71) 



where $ in Eq. ([TT]) means that if /3^p > /3l^^^ there 
is no stable boundary with the P region. A schematic 
representation of this case is given in Fig. 3. 

Case (4): If 2 < do < oo, /9c2 satisfies Eq. ((66)) (or, in 
terms of temperature Eq. (I67p for Tc2); whereas for (3ci 
we have either 



(3. 



cl 



(SG» 



(72) 




FIG. 1: Phase diagram for the case (1): Jo > and do < 2 
or do = cxD in a broad sense. 



if 3pcl ^ Pel > P, 



(SG) 



if 



>(3, 



(SG) 



(73) 



where in Eq. (I72|) we have introduced the symbol SG> to 
indicate that in general the stable P-SG surface is above 
(or below in terms of temperatures) the surface coming 

from the solution with label SG: /^f > /jf . Notice 
that, similarly to the case (3), we cannot exclude that 
Pel in Eq. (|73p be a non single- value function of X, as 

well as /3c in Eq. (|72p . A schematic representation 
of this case is given in Fig. 4. 

If more than two critical temperatures are present, the 
above scheme generalizes straightforwardly. 

Keeping our definition for the introduced symbol "SG 
and/or F" , we stress that: in all the fours cases the phases 
F and "SG and/or F" are exactly localized; in the cases 
(1) and (3) the phases P and SG are exactly localized; 
in the cases (2) and (4) the SG phase is always limited 
below (in terms of temperatures) by the unstable P-SG 
surface coming from the solution with label SG (indicated 
as P-SG unst in Figs. 2 and 4). Finally, we stress that 
- under the hypothesis that Co has only loops of even 
length - the stable P regions correspond always to the 
solution with label F. 

For 2 < do < oo, from the second hue of Eqs. (p^ 
and (|7^ and from Eq. ([T^]) . we see that the method is 
not able to give the complete information about the P- 
SG boundary since we have only inequalities, not equal- 
ities. Furthermore, in these regions of the phase dia- 
gram the critical temperature in general may be a non 



single- value function of X. On the other hand, we have 
the important information that in these equations the in- 
equalities between Tc and Tc (the solution with label 
SG) are always strict. As a consequence, we see that, 
when 2 < do < oo, in these regions the SG "magnetiza- 
tion" m^^^' will always have a finite jump discontinuity 
in crossing the surface given by Tc- In other words, along 
such a branch of the critical surface corresponding to the 
second line of Eqs. ^ and ^ and Eq. we have 

a first order phase transition, independently of the fact 
that the phase transition corresponding to the Ti^^^ sur- 
face is second-order, and independently on the sign of Jq. 



E. Generalizations 

The generalization to the cases in which the unper- 
turbed model has an Hamiltonian Hq involving couplings 
depending on the bond 6 € Fq is straightforwardly. In 
this case we have just to substitute everywhere in the 
formulae HHD-dSll), J^^^ with the set {d^^^}. However, 
the critical behavior will be in general different and more 
complicated than that depicted in the Subsections IIIBl 
and IIIB2. In particular, even in the case in which all the 
couplings Jgj are positive, we cannot assume that the 
Landau coefficient b'^^^ be positive so that, even in such 
a case, first-order phase transitions are in principle pos- 
sible, as has been seen via MC simulations in undirected 
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FIG. 2: Phase diagram for the case (2): Jo > and 2 < 
do < oo. 



FIG. 4: Phase diagram for the case (4); Jo < and 2 < 
do < oo. 




FIG. 3; Phase diagram for the case (3): Jo < and do < 2 
or do = oo in a broad sense. 



small- world models [l9[- 

As anticipated, our method can be generalized also to 
study possible antiferromagnetic phase transitions in the 
random model. There can be two kind of sources of anti- 
ferromagnetism: one due to a negative coupling Jq in the 
unperturbed model, the other due to random shortcuts 
Jij having a measure d/i with a negative average. 

In the first case, if for example the sublattice Cq is 
bipartite into two sublattices /Cg"^ and C^q'\ the unper- 
turbed model will have an antiferromagnetism described 



by two fields TOq"'' and TOq"''. Correspondingly, in the 
random model we will have to analyze two effective fields 
to'^''-' and to'^'') which will satisfy a set of two coupled self- 
consistent equations similar to Eqs. (|13p and involving 
the knowledge of mg"'' and mg^' . More in general, we can 
introduce the site-dependent solution toqi to find corre- 
spondingly in a set of coupled equations (at most N) , the 
effective fields of the random model. 

In the second case, following (2^ we consider a lat- 
tice Cq which is composed of, say, p sublattices C(^\ 
V = 1, . . . ,p. Then, we build up the random model with 
the rule that any shortcut may connect only sites belong- 
ing to two different sublattices. Hence, as already done 
in 13] for the generalized SK model, we introduce p effec- 
tive fields m^"^) which satisfy a system of p self-consistent 
equations involving the p fields mg"'' and calculated in 
the p external fields J^^^m*^"' (note that here the symbol 
F stresses only the fact that the effective coupling must 
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be calculated through Eq. (fT4)) ). 



IV. MAPPING TO NON RANDOM MODELS 

In Sec. V we will derive the main result presented in 
Sec. III. To this aim in the next subsection wc will recall 
the general mapping between a random model, built up 
over a given graph, and a non random one built up over 
the same graph, whereas in the following second subsec- 
tion we will generalize this mapping to random models 
built up over random graphs. We point out that the 
mapping does not consist in a sort of annealed approxi- 
mation. 



A. Random Models defined on Quenched Graphs 

Let us consider the following random model. Given a 
graph g, which can be determined through the adjacency 
matrix for shortness also indicate by g = {gb\^ with gb = 
0, 1, 6 being a bond, let us indicate with Vg the set of the 
bonds b oi g and let us define over Tg the Hamiltonian 



H ({(Ti}; {Jb}) = - 51 Jb(Ti,<Tj, - 

beFg I 



(74) 



where J\, is the random coupling at the bond 6, and 
tJij , (Tij are the Ising variables at the end-points of h. The 
free energy F and the physics are defined as in Sec. II 
byEqs. @-®: 



-HF- J dV {{Jb}) log (Z {{Jb})) , (75) 
(0)T= / dV{{Jb}){0)\ 1 = 1,2 (76) 



where dV {{Jb}) is a product measure over all the possible 
bonds b given in terms of normalized measures dfib > 
(we are considering a general measure dfj,b allowing also 
for a possible dependence on the bonds) 

dV {{Jb}) = n "^t^b (■^b) , / dfib {Jb) = 1, (77) 
berf„ii 

where Ffun stands for the set of bonds of the fully con- 
nected graph. As in Sec. II, we will indicate a generic 
correlation function, connected or not, by C with under- 
stood indices ii,...,ik all different, see Eqs. PT!]) and 

In the following, given an arbitrary vertex i of g, we 
will consider as first neighbors j of i only those vertices 
for which J d^j.ij{Jij)Jij or J d^ij{Jij)Jf ^ are at least 
0{l/N). Note that we can always neglect couplings hav- 
ing lower averages. We will indicate with D{Tg) the av- 
erage number of first neighbors of the graph g. For a 
d-dimensional lattice, D{rg) = 2d— 1, for a Bethe lattice 
of coordination number q, D{Tg) = q — 1, and for long 



range models, D{Tg) cx N. We will exploit in particular 
the fact that D{rco U Tf^n) = i:'(rfuii) oc N. 

Given a random model defined trough Eqs. (|74ll77p . 
we define, on the same set of bonds Tg , its related Ising 
model trough the following Ising Hamiltonian 

Hi ({^J; {Ji'^}) ^ - E Jb^^^^,3^ - E ^^'^^ (78) 



where the Ising couplings J^^^^ have non random values 
such that V 6, 6' e Tg 



jj,P = J^^' if d^b' = d^ib, 



jf'^ ^ if 



f dfib{Jb)Jb = O (jj) , 
J dfib{Jb)J^ ^ O {^) ■ 



(79) 
(80) 



In the following a suffix / over quantities such as Hj, Fj, 
fi, gi, etc. . . , or J^^\ /3c^\ etc. . . , will be referred to the 
related Ising system with Hamiltonian ([75)) . 

We can always split the free energy of the random sys- 
tem with N spins as follows 

- PF = J2 [ df^b log [cosh(/?Jfc)] + 
b •' 

^log[2cosh(/3/i,)] +0, (81) 



(j) being the high temperature part of the free energy. Let 
Lp be the density of cj) in the thermodynamic limit 



if = lim (p/N. 



N 



(82) 



Let us indicate with ipj the high temperature part of 
the free energy density of the related Ising model defined 
through Eqs. (|78ll80p . As is known, (pi can be expressed 
in terms of the quantities Zb = tanh(/3Jj'^') and Zi = 
tanh(/3/ii), i.e., the parameters of the high temperature 
expansion: 



ipi = ipi ({tanh(/3J,(^')};{tanh(/3;i,)} 



(83) 



The related Ising model is defined by a set of, typically 
few, independent couplings {J^^^}, trough Eqs. (|79ll80p 
and, for hi = 0, i = 1,...,N, its critical surface will 
be determined by the solutions of an equation, possibly 

vectorial, Gj ({tanh(/3jy^)}) = 0. 

In [l3| we have proved the following mapping. 

Let Z?!^^^ and f3c^^^^^ be respectively solutions of the 
two equations 

G/ (^{y"rf/ifctanh2(/3(SG)j,)}^ ^ 0, (84) 
Gi (^{ J d/i6tanh(/?f /AF) j^^ij ^ o_ (85^ 
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Asymptotically, at sufficiently high dimensions D{rg), 
the critical inverse temperature of the spin glass model 
/3c is given by 



/3. = min{/3f«),/?™} 



(86) 



and in the paramagnetic phase for D{rg) > 2 the follow- 
ing mapping holds 













C 



= o 



( 1 



(87) 



= \'^^^^ = y"^^ ( { / '^A*fotanh'(/3J;,)} 



and 



Ce// = jc(^) = jCi [{ I dA*btanh'(/3,A)}j , (89) 



where 



2, if ^idJdfibtanh'iPJb)}) 
2\ipi ({/ d^f,tanh(/3Jb)}) |, 
1, if ifi ({/ d^fetanh2(/3Jb)}) 
2\ipi ({/ d^ibtar^h{|3Jb)}) |, 



> 



< 



(90) 



and E=F/AF or SG, for l—l or 2, respectively. 



In the hmit D{Tg) 



and /i," 



0, 



1,...,7V, 



Eqs. (|84ll90p . give the exact free energy and correla- 
tion functions in the paramagnetic phase (P); the ex- 

act critical paramagnetic-spin glass (P-SG), /3c , and 

paramagnetic- F/AF (P-F/AF), /3c , surfaces, whose 
reciprocal stability depends on which of the two ones has 
higher temperature. In the case of a measure dfi not 
depending on the bond b, the suffix F and AF stand 
for ferromagnetic and antiferromagnetic, respectively. In 
the general case, such a distinction is possible only in 
the positive and negative sectors in the space of the pa- 
rameters of the probability distribution, {J d^i,Jb > 0} 
and {J d^bJb < 0}, respectively, whereas, for the other 
sectors, we use the symbol F/AF only to stress that the 
transition is not P-SG. 

It is not difficult to see that, when the measure dfi 
does not depend on the specific bond b, i.e., if dfXb = 
d^ib' V6, b' eTg, in the P region Eqs^^ (|84j|90p lead to the 
following exact limit for cp and C [15| 



lim (f 



hm C 

z5(rg)^oo 



0, for /3</3c, (91) 



therefore, the basic role of Eqs. (|87ll90p . is to show how, 
in the limit D(Tg) — > oo, and C approach zero and 
which are their singularities. In particular this proves 
that for all the (random) infinite dimensional models and 
any disorder non bond-dependent, the critical exponent 
a' for the specific heat has the mean-field classical value, 
a' = 0, and that the correlation functions (with differ- 
ent indices) above the critical temperature are exactly 



zero. We point out however that, when the measure d/if, 
depends explicitly on the bond b, Eq. (PT|) in general 
does not hold [23|. In fact, when the measure dfib is 
bond-dependent, the symmetry expressed by Eq. (j9ip is 
broken since the bonds are no longer equivalent. As we 
will see in the next section, in small-world models with 
an underlying lattice Cq having do < 2, even if Eq. (pij) 
may still holds for tp, the symmetry is broken for C since 
the direction(s) of the axis(es) of Cq is(are) now favored 
direction(s). Yet, if 2 < dp < oo, the symmetry (|91l) for 
if is broken as well. 

The analytic continuation of Eqs. (|87ll90p to /3 > /3c 
and/or for h ^ provide certain estimations which are 
expected to be qualitatively good. In general such es- 
timations are not exact, and this is particularly evident 
for the free energy density of the SG phase. However, 
the analytic continuation for the other quantities gives a 
good qualitative result and provide the exact critical be- 
havior (in the sense of the critical indices) and the exact 
percolation threshold. 

For models defined over graphs satisfying a weak def- 
inition of infinite dimensionality, as happens on a Bethe 
lattice with coordination number g > 2, a more general 
mapping has been established [l^. In this case, all the 
above equations - along the critical surface (at least) - 
still hold exactly in the thermodynamic limit, where we 
can set effectively D{rg) = oo. However, for the aims of 
this paper we do not need here to consider this general- 
ization of the mapping. 

We have yet to make an important comment about 
Eqs. ([23]), ((24|) and ([26]) . concerning the evaluation of 
a correlation function in the SG phase here for a ran- 
dom system with Jq = (for the moment being). In fact 
Eq. (|89p. for both a normal and a quadratic correlation 
function C'^^-' or C'^^\ has a factor 1/2 not entering in 
the physical Eqs. ^ and The difference 

is just due to an artefact of the mapping that separates 
the Gibbs state into two pure states |21j not only in the 
F case, but also in the SG case. In fact, let us con- 
sider the correlation functions of order fc = 1, that is, 
C(i' = and C(2) = QEA = We see that, for 

C(i), Eq. dlSl) in the SG phase gives C^^) = m'-^'^'^ /2. 
On the other hand, for any non zero solution to'^^^ of 
the self-consistent Eq. (fT3|) . there exists another solu- 
tion —m,^^'^\ and both the solutions have 1/2 probabil- 
ity to be realized in the random model. Since the SG 
phase is expected to be the phase characterized by having 
QEA 7^ and (ai) — 0, we see that if we introduce both 
the solutions ■m^^'^'> and — m^^*^-*, we get (cti) = in the 
SG phase. Similarly, for C^^), Eq. 1^ in the SG phase 
gives C'^' = (m(^*^')^/2, which at zero temperature gives 
1/2, whereas a completely frozen state with qea = 1 is 
expected. Again, we recover the expected physical qea 
by using both the solutions m*^^^^ and —m^^'^\ Repeat- 
ing a similar argument for any correlation function of 
order fc, and recalling that for k even (odd) the correla- 
tion function is an even (odd) function of the external 
magnetic field h, we arrive at Eqs. ([M)) and 
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B. Random Models defined on Unconstrained 
Random Graphs 

Let us consider now more general random models in 
which the source of the randomness comes from both the 
randomness of the couplings and the randomness of the 
graph. Given an ensemble of graphs g € G distributed 
with some distribution P{g), let us define 

Hgi{cr^}■,{Jb}) = - ^ JbfTi.CTj^ - /l ^ CTi 

beFg i 

The free energy F and the physics are now given by 

-f^F^Y. P^a) [ dV ({ JJ) log {Zg {{.h})) , (93) 
gee 

and similarly for (O). Here Zg{{J},}) is the partition 
function of the quenched system onto the graph realiza- 
tion g with bonds in Tg 



Zgi{Jb}) = J2- 



-/3ffg({<T,};{J.}) 



(94) 



and dp ({ J&}) is again a product measure over all the 
possible bonds b given as defined in Eq. ([77| . Note 
that the bond-variables {gb} are independent from the 
coupling- variables {Jt}- 

For unconstrained random graphs, or for random 
graphs having a number of constrains that grows suffi- 
ciently slowly with N, the probability P{g), for large N, 
factorizes as 



Pi9)= n Pbiab). 



(95) 



In such a case we can exploit the mapping we have pre- 
viously seen for models over quenched graphs as follows. 
Let us define the effective coupling Jf,: 



Jb = Jb ■ gb, 



(96) 



correspondingly, since the random variables Jb and gb are 
independent we have 



(97) 



dflb{Jb) = d^b{Jb) -Pbigb), 
with the sum rule 

dh{Jb)f{Jb\gb)= Y Pbi9b) / diJ.b{Jb)f{Jb;gbpS) 

36=0,1 

As a consequence, if we define the following global mea- 
sure 

dV ({ Jb}) = P{g) ■ dV {{Jb}) = n dh{Jb), (99) 



we see that the mapping of the previous section can be 
applied as we had a single effective graph Tp given by 

Lp = {& e Lfuii : Pb{gb - 1) ^ 0}, (100) 

in fact we have 



/3F = I dV [{Jb}) log (Zp ({.lb} 



(101) 



where Zp is the partition function of the model with 
Hamiltonian Hp given by 

Hp ({cT.}; {Jb}) =-J2 •^"'^'^'^J^ -hJ2''^- (102) 

6er„ i 



V. DERIVATION OF THE SELF-CONSISTENT 
EQUATIONS 

By using the above results, we are now able to derive 
easily Eqs. ()13ll2ip . Sometimes to indicate a bond b we 
will use the symbol (i,j), or more shortly ij. 

It is convenient to look formally at the coupling Jq also 
as a random coupling with distribution 

d/.o(4)/d4 = ,5(4 - Jo). (103) 

Let us rewrite explicitly the Hamiltonian ^ as follows 

He = - ^ {Jq + Ci-jJi-j) aiffj 
(ij)6ro 

- Y CijJijUiaj - h'Y'^i, (104) 

i<j, (*,j)^ro i 

and let us introduce the random variables J^, g'^^ and J^, 
where 



T/ _ ] Ja + CbJb, & e Lq, 
■""-^ Jb, bi Fo, 



9b 



Cb, 



and 



•Jb = 'Jb ' 9b- 



(105) 



(106) 



(107) 



Taking into account that the random variable Jq -I- Cij Jij , 
up to terms 0(1/A^), is distributed according to J^o(Jo), 
the independent random variables J^ and 5^ have distri- 
butions d/i'f, and p'j^ respectively given by 



and 



dfi'M) 



P'b{9b) = 



dfJ-oiJb), be To, 
d/i(J^), b^To, 



pig'b)^ 



be To, 



(108) 



(109) 
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where the measures dfj, and p are those of the model 
introduced in Sec. II. As a consequence, Eq. (|104p can 
be cast in the form of Eq. (jl02p with the measure 



(110) 



FinaUy, since Pb{gb) 7^ for any b e Ftun, we have also 



full, 



(111) 



and due to the fact that £'(rfuii) oc N, in the thermody- 
namic limit the mapping becomes exact. 

According to Eqs. ()78ll80p . the related Ising model 
of our small-world model has the following Hamiltonian 
with two free couplings: Jq^\ for Fq, and J^^\ for Ffun 



Hi - -Jg 



(112) 



After solving this Ising (/) model the mapping allows 
us to come back to the random model by performing 
simultaneously for any b G Ffun the reverse substitutions 

tanh ^ f dfL',{4) (4) tanh' (pj',) , (113) 



where ^ = 1, 2 for S = F or SG solution, respectively. 
Since the couplings Jq^' and J*^^) are arbitrary, we find it 
convenient to rcnormalize J*^'^^ as j'^^^ /N and at the end 
of the calculation to put again J(^) instead of jt-^'/TV. 
Note that for the mapping nothing changes if we do not 
make this substitution; the choice to use j'^' /N instead 
of j'^' is merely due to a formal convenience, since in 
this way the calculations are presented in a more stan- 
dard and physically understandable form. In fact, ac- 
cording to Eqs. (jllOp and (|113p what matters after solv- 
ing the related Ising model with J^^^ /N instead of J^^^ 
is that, once for S=F and once for S=SG, we perform 
- simultaneously in the two couplings - the following re- 
verse mapping transformations (I — 1,2 for E = F or SG, 
respectively) : 



tanh (/3J(^V^) J djl (jy) tanh' (/3Jy) (114) 
for ^ To, and 



tanh 



Jij I tanh 



, (115) 



for e Fq. 

Explicitly, by applying Eqs. piOp and Eq. ^ the 
transformations (|114p and (|115p become, respectively 



and 



PJq > i3Jq ^ , 



(116) 
(117) 



where we have made use of the definitions (fT4|) - ((T7)) in- 
troduced in Sec. III. 

Let us now solve the related Ising model. We have to 
evaluate the following partition function 

Zi = g/3"'y'E(i,j)ero '^>'^3+/34r-2:i^j<^>'73+/3/iE>f^,«218) 

In the following we will suppose that J*^^) (and then J^^^) 
is positive. The derivation for J^^^ (and then J'-^-') neg- 
ative differs from the other derivation just for a rotation 
of 7r/2 in the complex m-plane, and leads to the same 
result one can obtain by analytically continue the equa- 
tions derived for J^^^ > to the region J*^^) < 0. 

By using the Gaussian transformation we can rewrite 
Zj as 



Zi = cn^^ e 



(>,j)ero ' 



dm e-ij'-''rn^N+p{j^'^m+h) (119) 

where cat is a normalization constant 



Cat 



2tt 



(120) 



and, in the exponent of Eq. (|119p . we have again ne- 
glected terms of order 0(1). For finite N we can exchange 
the integral and the sum over the cr's. By using the def- 
inition of the unperturbed model with Hamiltonian Hq, 
Eq. JT]), whose free energy density, for given PJq and 
is indicated with fodSJo^/Sh), we arrive at 



Zl = CN 



dm e-^^("), 



where we have introduced the function 



^ ^ 2 



(121) 



(122) 



By using dph PfoiPJo,l3h) = -mo{f3Jo,l3h), and 
d/sh mo{pjo,ph) = xo{(3Jo,(3h) we get 

i'(m)=/3j(^) m~ma(^(34^\(3j'-^^m + f3hj ,(123) 
L"{m) = (3J^^'> 1 - /3j(-^^xo (^(34'\(3J^''>m + Jl24) 



If the integral in Eq. (|12ip converges for any N, by per- 
forming saddle point integration we see that the saddle 
point m'^P is solution of the equation 

m'P = mo (/3J^^\/3j(^)m"P 4-/3/i) , (125) 

so that, if the stability condition 

1 - 13J^^'>XQ (/3J^^\/3j(-^)m^P + Ph^ > 0, (126) 
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is satisfied, in the thermodynamic hmit we arrive at the 
foUowing expression for the free energy density // of the 
related Ising model 



In the region of temperatures where Eq. (|129p is vio- 
lated, Eq. (|125p gives two symmetrical stable solutions 
±m/ ^ 0. From Eq. (|130p we see also that the case 
l3c'^^ — (3\^Q is impossible unless be J^^^ = 0, since the 



J(^)to^ + /3/o (/3 J^^^ , /3 J(^)to + [3h^ (127) susceptibility xoiPJo'' : 0) must diverge at /3^o^ . We have 



Similarly, any correlation function C/ of the related Ising 
model is given in terms of the correlation function Cq of 
the unperturbed model by the following relation 



(3h 



(128) 



Of course, the saddle point solution m^P represents 
the magnetization of the related Ising model, as can be 
checked directly by deriving Eq. (|127p with respect to 
I3h and by using Eq. 

If the saddle point equation p25p has more stable solu- 
tions, the "true" free energy and the "true" observable of 
the related Ising model will be given by Eqs. (|127p and 
P28p . respectively, calculated at the saddle point solu- 
tion which minimizes Eq. (|127p itself and that we will 
indicate with m/. 

Let us call [3^!^ the inverse critical temperature of the 
unperturbed model with coupling J^^^ and zero external 
field, possibly with /3^g^ = cxd if no phase transition exists. 
As stressed in Sec. IIIB, for the unperturbed model we 
use the expression "critical temperature" for any temper- 
ature where the magnetization mg at zero external field 
passes from to a non zero value, continuously or not. 
Note that, as a consequence, if Jg'^'' < 0, we have for- 
mally p'^J^ = oo, independently from the fact that some 
antiferromagnetic order may be not zero. 

Let us start to make the obvious observation that a 
necessary condition for the related Ising model to have a 
phase transition at ft. = and for a finite temperature, 
is the existence of some paramagnetic region P/ where 
TO/ = 0. We see from the saddle point equation l|125p 
that, for ft, = 0, a necessary condition for to/ = to 
be a solution is that be /9 < (3^^^ for any in P/, from 
which we get also (3^^^ < /3^q'' . In a few lines we will see 

however that the inequality must be strict if /3^q' is finite, 
which in particular excludes the case Jq < (for which 
the inequality to be proved is trivial). 

Let us suppose for the moment that be /3c '^^ < /3^q\ 

For (3 < P^Q^ and ft = 0, the saddle point equation (|125p 
has always the trivial solution mj = which, according 
to Eq. (|126p . is also a stable solution if 

1-PJ^'^Xo{p4'\o) >0- (129) 
The solution nij — starts to be unstable when 

1-/3J'-'^Xo[p4'\o)^0- (130) 

Eq. (|130p , together with the constrain (3^^^ < (3^^^ , gives 
the critical temperature of the related Ising model I3i^\ 



therefore proved that (3^^^ < (3^^^ . Note that for J^^^ > 
and (3 < ^^0^ Eq. is violated only for /3 > (3i^\ 

whereas for Jq^'' < Eq. (|129p in general may be vio- 
lated also in finite regions of the (3 axis. 

The critical behavior of the related Ising model can be 
studied by expanding Eq. (|125p for small fields. How- 
ever, we find it more convenient to expand L{m) in series 
around to = since in this way everything can be cast 
in the standard formalism of the Landau theory of phase 
transitions. From Eq. (|122p . taking into account that 
the function xo {(3Jo,l3h) is an even function of /3ft, we 
have the following general expression valid for any to, (3 
and small h 



(m) = /3/o (/3 J^'^ , O) - TOO (/? J^^ , O) /3ft {mllil) 



L 



where we have introduced the Landau free energy density 
'il'{m) given by 

lb (m) — —ara^ H — brn^ H — cto^ 
2 4 6 

-TO/3ft + A (/3/o) (/3jy\/3j(^)m) , (132) 



where 



/3j(^)xo(/5^r,0)]/3J(^\ (133) 



9(/3ft)2 ^" 



d{l3h)' 



■Xo 



3! 



J)) 



5! 



, (134) 



-, (135) 



ft = TOO (/3J^'\o) J^''>+xo {p4'\o) /3(^)j(^)/3ft(136) 

finally, the last term A (/3/o) {j34'\ (3J^^^m^ is defined 

implicitly to render Eqs. (|13ip and (|132p exact, but 
terms 0{h^) and Oim^h); explicitly 



A 



(/3/o)(/3Jr,/3J('' 



E 

k=4 



Q2k-2 

2k^ ^0 



diPh) 



, 2k 



lih = Q 



{2k)\ 



4137) 



Finally, to come back to the original random model, 
we have just to perform the reversed mapping transfor- 
mations (fTTB)) and (fTTTl) in Eqs. P^ - (fTT7)) . As a result 
we get immediately Eqs. (II31)-(|1I]), but Eq. (dH). 
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VI. DERIVATION OF EQ. ([19]) AND EQS. 



Concerning Eq. p9)) for the full expression of the free 
energy density, it can be obtained by using Eqs. ([5T|) . 
([5^ . ((55)) and (HOI)- Here (^/ is the high temperature 
part of the free energy density of the related Ising model 
we have just solved: 



f3fi = Um ^ log cosh(/34 



N -1 



(«j)Gro 



X log 



cosh 



(/3J(^V^)] +log[2cosh(/3/i)] +<y97(138) 



where we have taken into account the fact that our re- 
lated Ising model has jFol connections with coupling J^P 
and N{N — l)/2 connections with the coupling j'^'/A^. 
By using Eq. ()127l) calculated in nii and Eq. (I138|) . for 
large TV we get 

1 



lim — log cosh(/3Jn^h 



(i,i)Gro 



log[2cosh(/3;i)]+0( - 



(139) 



Therefore, on using Eq. ([55)) . for the non trivial part 
93^^) of the random system, up to corrections 0(1/7V), 
we arrive at 



where 



ipoipJo,l3h)^-l3fo(^l3jP,ph 



1 



lim ,^ 

Af^oo TV 



^ log[cosh(/3jf^; 
log [2cosh(/3/i)] , 



(143) 



is the high temperature part of the free energy density 
of the unperturbed model with coupling Jq and exter- 
nal field h. There are some important properties for the 
function (po(/3Jo,0): it is a monotonic increasing function 
of pjo; if the lattice Co has only loops of even length, 
ipo{PJo,0) is an even function of /3Jo; furthermore, if 
do < 2 and the coupling-r ang e is finite, or ii do = oo 
at least in a wide sense [l5|, in the thermodynamic 
limit we have (poiPJojO) = 0; if instead 2 < do < oo, 
fo{PJo,(^) 7^ 0. We see here therefore what anticipated 
in Sec. IVA: when Jq ^ 0, the symmetry among the ran- 
dom couplings is broken and for do sufficiently high this 
reflects in a non zero (p^^^^ also in the P region. 

Next we prove Eqs. ([64]) -([73 ]) . To this aim we have 
to calculate Eq. (|142p at the leading solution m*^^) and 
to compare and ip'^^'^\ Note that the term in the 
square parenthesis of Eq. (|142[) is non negative since 
m*^^^ is the absolute minimum of L(^) . We recall that 
for critical temperature we mean here any temperature 
lying on the boundary P-F or P-SG, so that 
for any (3 in the P region. 



Jo > 



^(S) = (ni^^)y - log [2 cosH(3h)] 

— lim — log cosh(/3Jn^'') 

(i,i)ero 

-/3/o(/?jf\/3J^''^m(^)+/3/i) (140) 
In terms of the function L(^^(m) Eq. (jl40p reads as 



hm 1 J2 cosh(/3jf)) 
log [2cosh(/3/i)] . 



(141) 



By using Eqs. ([51]), piT|) . ([SS]) and with ^ 1 or 2 
for Ti—F or E^SG, respectively, we get Eq. (jl9p . 

For h = Eq. (|14ip can conveniently be rewritten also 

as 



(0) - l(^) (m(^)) 



(142) 



If Jo > 0, for both the solution with label F and SG, 
we have only one second order phase transition so that 
7Ti(i^) = and TO^^^) = 0, respectively, are the stable 
and leading solutions even on the boundary with the P 
region. 

Let us suppose /3c^'' < /Jc^*^'. Let be ifoi-, 0) ^ 0. From 



Eq. p42p and by using Jq^^ > Jg^"^'' , we see that 



.(SG) 



> 



(F) 



^0 ^/3r4' 

^(^«)i,<.,=^o(/3rjr\o) 



(144) 



Finally, by using this result and the general rule given by 
Eqs. ((88)) and ([90]) . we see (and with a stronger reason, 
due to the factor 1/2 appearing in these equations for 
the SG solution) that the stable phase transition is the 



P-F one: f3c = (3^ . Similarly, by using Eq. ITi^ for 

(F) (SG) 

(3c < f3 < (3c , we see that even for any in the 
interval (/3c^'' , /Si^*^'' ) the stable solution is that with label 
F. This last observation makes also clear that if fo{', 0) — 
wc reach the same conclusion: F is the stable phase in 
all the region /Sc^' < (3 < (3^^^ and in particular this 
implies also that the stable phase transition is the P-F 

one: (3c = (3c ■ 
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Let us suppose Pc^^ > /?c^^^ If v9o(-,0) ^ 0, we arrive 



at 



(SG) 



(145) 



Finally, by using this result and the general rule given 
by Eqs. and we see (and with a stronger rea- 
son) that the stable phase on the boundary is that pre- 
dicted by the F solution which has zero magnetization at 
^(SG) rpj^jg does not imply that /3, 

If instead (fo{-,0) 



CP) 

(3c , but only that 



0, by using Eq. 
(fT42l) for Pf"^' < (3 < P^''', we see that for any /3 in the 
interval {pi^'^\pi^^) the stable solution is SG and then, 



(SG) 



in particular, the stable boundary is P-SG: /3c = Pi 



B. Jo < 

If Jo < 0, for the solution with label F, we may have 
both first and second order phase transitions. In the first 
case we cannot in general assume that is the stable 
and leading solution on the boundary with the P region: 



(F) ^ in general. As a consequence, for a first or- 



der transition the term in square parenthesis of Eq. (jl42p 
may be non zero even on the critical surface. Further- 
more, as Jo < 0, for the solution F we have at least two 
critical temperatures that we order as P^^ > P^2 ■ How- 
ever, despite of these complications, if we assume that Cq 
has only loops of even length, (po{-^ 0) turns out to be an 



even function and, due to the inequality |Jo^^| > 
almost nothing changes in the arguments we have used 
in the previous case Jo > 0. 



.(SG) 



Let us consider first the surfaces P^2 ^'^d P, 



,(SG) 



Inde- 



pendently of the kind of phase transition, first or sec- 
ond order, we arrive again at Eqs. (|144p and (|145p . 

for /3^2^ < Pc^'^^ and /3^2^ > pi^'^\ respectively, with 
the same prescription for the cases ipo{-,0) ^ 0, or 
<^o(-,0) = 0. 

Let us now consider the surfaces /3^P and Pc^'^\ If 



/3^P < Pc°^' and (po{-,0) ^ 0, for any P in the interval 
[P^^\pi^^^] we have. 



,(SG) 



{p4 



(F) 



>^(s«)|, = ^o {pjr\0) 



(146) 



so that the interval [P^^^ , pi^'^^] is a stable P region corre- 
sponding to the solution with label F. Similarly, we arrive 
at the same conclusion if (po(-,0) = 0. However, the in- 
terval of temperatures where the P region [p'^\ pi^'^^] 
is stable can be larger when do > 2. In fact (exactly 
as we have seen for Jo > 0) in this case the P-SG stable 
boundary may stay at lower temperatures. Finally, let us 



consider the case /J^i > Pc ■ If </?o(-, 0) = 0, by using 
Eq. (|142p we see that for any P > /^c^*^' we have that 
the stable solution corresponds to the SG one, so that 
there is no stable boundary with the P region. If instead 
ipo{-,0) ^ 0, due to the fact that (p^^^^j^ and if^^^p grow 
in a different way with P, we are not able to make an ex- 
act comparison, and it is possible that the P-F boundary 
becomes stable starting from some Pd with Pd > PcV ■ 
In general, as in the case Jq < 0, we could have one (or 
even more) sectors where the P region corresponding to 
the solution with label F is stable. 



VII. CONCLUSIONS 

In this paper we have presented a novel and general 
method to face analytically random Ising models defined 
over small- world networks. The key point of our method 
lies on the fact that, at least in the P region, any such 
a model can be exactly mapped on a suitable fully con- 
nected model, whose resolvability is in general non trivial 
for do > 1, but still as feasible as a non random model. 
The main result stands then in deriving a general self- 
consistent equation, Eq. (fO|) , able to describe effectively 
the model and based on the knowledge of mo{pJo, ph), 
the magnetization of the unperturbed model in the pres- 
ence of an arbitrary external field. 

The physical interpretation of the these results is 
straightforward: from Eqs. ([0)) we see that, concerning 
the magnetization m,^^\ the effect of adding long range 
Poisson distributed bonds implies that the system - now 
perturbed - besides the coupling Joi feels also an effective 
external field J*^^) shrunk by m'^^ itself; whereas, con- 
cerning m'^^'^\ the effect is that the system now feels a 
slightly modified effective coupling /JJg^^-* and an effec- 
tive external field J^^'^^ shrunk by to^^*^^ itself. 

We are therefore in the presence of an effective field 
theory which, as opposed to a simpler mean-field theory, 
describes m^^^ and m'^^^ in terms of, not only an effec- 
tive external field, but also in terms of the non trivial 
function mo{pJo, Ph) which, in turn, takes into account 
the correlations due to the non zero short-range coupling 
Jq or Jq^'^^ felt by the unperturbed model. The combi- 
nation of these two effects gives rise to the typical be- 
havior of models defined over small-world networks: the 
presence of a non zero effective external field causes the 
existence of a mean-field phase transition also in low do 
dimension. However, the precise determination of both 
the critical surface and the correlation functions is deter- 
mined in a non trivial way by the unperturbed magneti- 
zation molPJo, ph). 

In this paper we have then applied the method to an- 
alyze, in the most general feasible situation, the critical 
behavior of generic models with Jo > and Jo < 0, 
showing that they give rise to two strictly different sce- 
nario of phase transitions. In the first case, we have a 
mean-field second-order phase transition but with a fi- 
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nite correlation length. Whereas, in the second case, we 
have multiple first and second order phase transitions. 
Furthermore, we have shown that the combination of 
the F and SG solution, results in a total of four pos- 
sible kind of phase diagrams according to the cases (1) 
(Jo > 0; do < 2, or do = oo), (2) (Jo > 0; 2 < do < oo), 
(3) (Jo < 0; do < 2, or do = oo), and (4) (Jq < 0; 2 < 
do < 00 ). One remarkable difference between the cases 
with do < 2, or do = oo and the cases with 2 < do < oo, 
is that in the latter we have also first-order P-SG transi- 
tions and, moreover, re-entrance phenomena are in prin- 
ciple possible even for Jo > 0. 

In a forthcoming paper (part II) we will apply the 
method to study in detail several solvable (by our 
method) models of interest. 

As discussed in Sec. HIE, the method can be readily 
generalized to study more complex phase diagrams that 
may emerge when one considers couplings Jq depending 
on the bond b, and to study small-world antifcrromag- 
netism. 

Models defined on complex small-world networks (2^ 
are interesting subject of future works. 
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APPENDIX A: GENERALIZATION TO NON 
HOMOGENEOUS EXTERNAL FIELD 

In this appendix we prove Eq. (|3T|) calculating the 
0{l/N) correction responsible for the divergence of the 
susceptibility of the random system at Tc- To this aim we 
firstly need to generalize our method to an arbitrary ex- 
ternal field. Let us consider again a fully connected model 
having - as done in Sec. V - long-range couplings J (for 
brevity we will here omit the label /) and short-range 
couplings Jo but now immersed in an arbitrary (non ho- 
mogeneous) external field {/i„}, where n — 1, . . . , N. Af- 
ter using the Gaussian transformation we have the fol- 
lowing partition function: 



and 



dm e-^^(") 



(Al) 



where we have introduced the function 

L{m) = ^ Jm^ + /3/o {f3jQ,{f3Jm + /3/i„}) , (A2) 

fo {PJo,{P^n}) being the free energy density of the un- 
perturbed model in the presence of an arbitrary external 
field {(3hn}. By using 

dfjh, /3/o(/3Jo, {PK}) = -mo,(/3Jo, {/3/i„}), (A3) 



df3h, moi{f3jQ,{f3hn}) = Xo-i,j{f3Ja, {/Shn}) 

= (o-io-j)o - (o-i)o(CTi)o (A4) 



we get 
L'(to) = (3J 



^ J2 mo^ (/3Jo, {f3Jm + phn}) 



(A5) 



L"{m) = PJ 



1-/3 J 



iV 



(A6) 



By performing the saddle point integration we see that 
the saddle point rrf^ is solution of the equation 

m^P = l^mo,(/3Jo,{/3Jm^P + /3/i„}) (A7) 

i 

hence, by using 

Xo {l3Jo,{pJm + I3K}) = 



N 

we see that if the stability condition 

1 - pjxo iPJoAPJm'P + Phn}) > 0, 



(A9) 



is satisfied, in the thermodynamic limit we arrive at the 
following expression for the free energy density / of the 
related Ising model immersed in an arbitrary external 
field 



^ Jm^ + /3/o (/3Jo, {pjm + (3K}) 



(AlO) 



On the other hand, by derivation with respect to Phi and 
by using Eq. (|A7[1 . it is immediate to verify that 



m, = (a,) mo^iPJo, {/3Jm^P + pK}), (All) 

and then also (from now on for brevity on we omit the 
symbol ''p) 



(A12) 



We want now to calculate the correlation functions. 
From Eq. (lAlip . by deriving with respect to Phj we 
have 



dr 



diph,) 



J2^o-aAPJo,{PJrn + Phn}) 



, ^ ^ dm ^ 



(A13) 
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which, by summing over the index i and using (|A12p . 
gives 



dr 



d{l3hj) 1 - pJxoiPJoAPJm + I3hn}) 



(A14) 



We can now insert Eq. (|A14p in the rhs of Eq. (|A13|) to 
get 



Xi] = XO;i,i(/?-/o,{/3>/"^ + /3/l„}) 

N 1 - /3Jxo(/3Jo, Wm + /3/i„}) ' 



(A15) 



where for brevity we have omitted the argument in Xo:i,j 
and xo;i,k, which is the same of xq appearing in the de- 
nominator. If we now come back to choice a uniform 



external field /i„ = h, n = 1, . . . ,N, we can use trans- 
lational invariance and for the related Ising model (fully 
connected) we obtain the following correlation function 



X^J 



(3 J [xo(/3Jo,/3Jm + /3fe)]' 
N 1- f3Jxoi(3Jo,(3Jm + (3h) 
Xo;i,3{0Jo,f3Jm + Ph). 



(A16) 



Finally, by performing the mapping substitutions (|116D 
and (fTTT]) we arrive at Eq. ([3T|) . 

Similarly, any correlation function C of the related 
Ising model will be given by a similar formula with the 
leading term Co plus a correction 0(1/A^) becoming im- 
portant only near Te- 
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